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ABSTRACT 



Selective Reenlistment Bonuses (SRBs) are offered to 
improve retention in designated military occupational 
specialties (MOSs) for specified years-of-service intervals 
(zones). The amount of the bonus is set by assigning an "SRB 
Multiplier" for each MOS and zone combination (cell). 

Determination of multipliers is modeled as a nonlinear 
knapsack problem which is then linearized to a generalized 
assignment problem. The objective is to minimize the sum 
over all cells of a weighted squared deviation from the 
reenlistment target in each cell. Lagrangian relaxation 
provides lower bounds and feasible solutions. The best 
feasible solution is improved using a greedy heuristic to 
apportion unexpended funds. 

A FORTRAN 77 computer program implements the procedure. 
Data for FY86 yields a 0-1 integer program with 4895 binary 
variables and 980 constraints. A solution within .01% of 
optimality is obtained on an IBM 3033AP in 1.7 seconds and 
on an IBM PC in about four minutes. 
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I. INTRODUCTION 



The United States Marine Corps (USMC) uses Selective 
Reenlistment Bonuses (SRBs) to increase retention in desig- 
nated military occupational specialties (MOSs) and years-of- 
service intervals (zones). Those MOSs and zones which will 
receive bonuses as well as the dollar amounts of the bonuses 
are determined by the assignment of "SRB Multipliers." In 
this study, the problem of assigning multipliers given a 
limited budget is modeled and solved as a mathematical 
program. In this chapter, the SRB program is described, 
terms are explained and the problem itself is defined. In 
addition, the solution approach is briefly discussed and the 
structure of this thesis is outlined. 

A. PROBLEM STATEMENT 

The Marine Corps SRB program was developed as a reen- 
listment incentive to increase retention in designated MOSs. 
In effect, an SRB is a sum of money offered to individuals 
with specific skills and years in service to encourage them 
to reenlist and thereby selectively increase manning levels. 
The bonus is applied at reenlistment points that fall 
between 21 months and 14 years of active duty service (ADS). 
This period is divided into three reenlistment zones. As 
defined in the applicable Department of Defense (DoD) direc- 
tive [Ref. 1], the period between 21 months and 6 years of 
ADS is designated as Zone A, the period between 6 and 10 
years as Zone B and the period between 10 and 14 years as 
Zone C. 

SRBs are not offered to everyone who reenlists, rather 
they are directed at selected MOS and zone combinations, 
called "cells" in this study. For each cell, there exists a 
reenlistment target which is defined as the desired number 
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of reenlistments in that cell. In many cells, enough Marines 
will reenlist without any added inducement so that the reen- 
listment targets will be automatically met or exceeded. In 
others, however, the absence of a bonus will result in 
insufficient reenlistments. Through the offer of bonuses to 
reenlist, it is possible to reach or approach the targets. 
To permit additional flexibility in awarding bonuses, the 
amount of an SRB can vary from a minimum of 0 to a maximum 
of $16,000. From experience, it is known that larger 
bonuses result in greater retention rates although the rate 
at which retention increases with the size of the bonus 
differs greatly between cells. The total amount of bonus 
money that can be awarded per fiscal year is limited by the 
SRB program's annual budget. Since an individual receives 
X5% c^f the bonus as a lump sum payment with the remainder 
apportioned in equal annual payments over the term of the 
reenlistment, the total SRB budget must consist of two 
parts. The first is money owed from bonuses from previous 
years and its apportionment is non-discretionary . The rele- 
vant budget here is the remainder which can be allocated for 
current year bonuses. 

Once the decision is made to consider the use of a non- 
zero SRB to improve manning in a particular cell, the actual 
monetary size of the bonus is based on increments of the 
eligible member's monthly base pay. The individual receives 
an amount equal to his monthly base pay multiplied by that 
cell's SRB multiplier and the number of years of additional 
obligated service, subject to a maximum of $16,000. The DoD 
directive stipulates that the multipliers may not exceed 6 
and may be fractional amounts. The Marine Corps' implementa- 
tion of this policy is somewhat more restrictive. It 
chooses not to use fractional multipliers and does not 
permit a multiplier greater than 5 for reenlistments in Zone 
A, 4 in Zone B and 3 in Zone C. Assigning a multiplier of 0 
is equivalent to no bonus . 
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Presently, SRB multipliers are assigned man ually in an 
iterative procedure. Multipliers are selected for cells 
based on the desired number of reenlistments for those 
cells. Once a particular combination of SRB multipliers has 
been selected, a microcomputer based spreadsheet package is 
used to determine if the proposed combination is feasible, 
i.e., within the allowed budget. If it is not, multipliers 
are changed and the calculation is repeated. 

At this time, there is no objective function with which 
alternative sets of proposed multipliers can be evaluated or 
compared. Instead, the officer responsible for assigning 
multipliers relies exclusively on his experience, judgement 
and various rules of thumb to judge the desirabililty of a 
particular set of multipliers. Even if a sensible objective 
function were defined, it is doubtful that an optimal solu- 
tion could be determined using a manual procedure. Of the 
approximately 1000 cells, on average, 450 are assigned non- 
zero SRB multipliers during any given period. It is very 
likely that the dimensionality of the problem would quickly 
overwhelm any ability to manually arrive at an optimal 
allocation . 

This study proposes an alternative procedure. It is 
automated due to the large amount of data. An objective 
function is defined and justified, thereby permitting the 
comparison of different sets of multipliers. Further, an 
optimization method built around this objective function and 
the budget constraint is described and tested. Using this, a 
nearly optimal set of SRB multipliers can be determined. 
This will result in a more effective allocation of the SRB 
program budget. 

B. PROBLEM SCOPE 

The goal of the SRB program is to reduce expected short- 
falls in the number of reenlistments in particular cells by 
increasing the retention rate in those cells through the 
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offer of a reenlistment bonus. The model that is developed 
in this study tries to achieve that goal by determining a 
set of SRB multipliers which, by a selected definition , can 
be considered optimal. This is a single period model; no 
attempt is made to forecast sets of multipliers for subseq- 
uent years based on expected retention resulting from the 
current year’s set of multipliers. It allocates that part 
of the budget not obligated for annual payments for bonuses 
from previous years and does not consider the effect of the 
25% of the bonus that is apportioned in future years. 
Furthermore, the model, through its use of expected values, 
is essentially deterministic. Other assumptions regarding 
the model are explicitly addressed in the succeeding 
chapter . 

C. MODEL AND SOLUTION APPROACH 

In the mathematical formulation that is proposed to 
model the assignment of SRB multipliers, the objective func- 
tion to be minimized is the sum of individual cell penal- 
ties. These penalties are nonlinear increasing functions in 
the deviation from the reenlistment target for the cell. 
Inclusion of the budget constraint, which is also nonlinear, 
results in a nonlinear knapsack problem. The nonlinearity, 
coupled with other complicating factors, prevents the 
problem from being efficiently solved using common knapsack 
algorithms. For descriptive and computational purposes, the 
problem is reformulated as a special generalized assignment 
problem. 

In this generalized assignment formulation, cells repre- 
sent the tasks that must be assigned. Unlike conventional 
generalized assignment problems, there is but a single agent 
who must accomplish all the tasks but who has alternative 
processes for each. These alternative processes correspond 
to the various multiplier levels permitted for that cell. 
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Because of the special structure of this problem, 
Lagrangian relaxation is particularly well suited as a solu- 
tion technique. Lagrangian relaxation has been successfully 
used in scheduling (e.g., [Ref. 2]), location 
(e.g., [Ref. 3]) and set-covering problems (e.g., [Ref. 4]) 
as well as generalized assignment problems [Ref. 9]. Of 
interest is that previous applications to generalized 
assignment problems have shared the characteristic of 
multiple tasks and multiple agents; a review of the litera- 
ture has not revealed an application with a formulation 
similar to the one developed in this study. 

An important advantage to the Lagrangian relaxation 
approach is that, in the process of establishing a lower 
bound on the optimal solution, feasible solutions are also 
obtained. A heuristic method is developed which improves 
the best feasible solution uncovered in the bounding proce- 
dure. The heuristically improved set of multipliers is 
accepted as the final solution provided its objective value 
is sufficiently close to the lower bound. 

D. THESIS OUTLINE 

This thesis develops and presents a method for deter- 
mining SRB multipliers in the USMC. In Chapter II, the 
problem is formulated first as a nonlinear knapsack problem. 
Then, to facilitate the description of the solution approach 
and for computational purposes, it is reformulated as a 
special generalized assignment problem. The solution meth- 
odology is presented in Chapter III. The approach uses the 
technique of Lagrangian relaxation combined with a heuristic 
procedure to provide a final solution. Details on the 
implementation of this process and computational results are 
presented in Chapter IV. Conclusions and recommendations 
are contained in Chapter V. Lastly, a listing of the source 
code, user instructions and samples from input and output 
files are included as appendices. 
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II. MODEL FORMULATION AND DESCRIPTION 



In this chapter, the process of determining multipliers 
for all cells is formulated first as a simple nonlinear 
integer program. Coefficients are defined and described and 
the rationale behind the selection of the objective function 
is explained. For explanatory and computational purposes, 
the problem is then converted to a special case of the 
generalized assignment problem. 

A. MODEL DEVELOPMENT 

The following formulation is developed to model the 
assignment of SRB multipliers: 

Indices : 

i = 1 , . . . ,m cells 



Data : 



B 



budget 



Functions : 




penalty for setting multiplier at x^^ 
in cell i 

cost of setting multiplier at x^ 
in cell i 

maximum permissible multiplier for 




cell i. n^ = 5,4,3 for Zone A,B,C 



respectively 



Decision Variables: 



multiplier for cell i 



Formulation: 



min Xc^(x^) 



(PI) 



s.t. ^ B 
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i = 1 , . . . , in 



0 ^ ^ 

integer 



In defining the objective function it is assumed that 
each cell's effect on the overall objective function is 

independent of all other cells. The objective function, 
then, is just the sum of individual cell penalties where a 
Ci(xi) should be 0 if the desired number of Marines are 

expected to reenlist given an SRB multiplier of x , and 

should rise as the expected deviation from the reenlistment 
target increases. An individual term in the budget 

constraint represents the expected resource requirements 
associated with selecting a particular multiplier value for 
that cell. For this formulation the following functional 
forms are selected: 



:i(xi) 



= (TiA^"lWi)x (eqn 2.1) 



ri(Xi) 



(R^(x^)E^) {min [ ( . 75P^S£X^) , 12000] } 



(eqn 2.2) 



where 

normalized cost of training a Marine in cell i 
A^ number of Marines in cell i 

exogenous cell weighting factor 
reenlistment target for cell i 
Ri(Xi) reenlistment response rate in cell i when 
multiplier x^ is offered 

E^ number of Marines in cell i eligible for reenlistment 

average monthly base pay of a Marine in cell i 
duration of reenlistment of a Marine in cell i 
Q relative weight of overages to shortages with 

respect to the reenlistment target 
(a-b)"^ positive part of (a-b) 
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The function Cj^(x^) represents a weighted squared devia- 
tion from the target number of reenlistments that would be 
expected in cell i if multiplier x^ were offered. It is the 
product of two sets of terms. The first term repre- 

sents a "weight" assigned to each cell. A weighting scheme 
is required as deviations in some cells are considered more 
critical than in others. The second set of terms 

{[(Di-Ri(Xi)Ei)^]2 + Q[(R. (xi)Ei-Di)^]2> 

is the squared deviation from the desired number of reen- 
listments in cell i but where Q may be used to weight the 
penalty incurred by overages (the second difference term) 
differently from the base penalty incurred by shortages (the 
first difference term). The Q, the "over/under" factor, has 
been set to 1 in this study but values somewhat less than 1 
might be appropriate. 

Squaring the deviation has the effect of penalizing 
large deviations much more severely than smaller ones. This 
is appealing because it evens out shortages and overages 
among cells. Most would agree that it is less detrimental 
to the force as a whole to "spread the shortages around" 
than to have a potentially debilitating shortfall in a few 
cells. Spreading overages around is desireable for logis- 
tical reasons . 

The deviation terms are calculated by squaring the 
difference between the desired number of reenlistments 
and the expected number of reenlistments R^^(Xj^)E^. The 
function Rj^(xj^) is discrete and defined only at integer 
values of x^^ on the interval [0,nj^]. Specifically, when the 
individuals in cell i are offered an SRB multiplier Xj^ , a 
fraction R^^(x^) of the E^ individuals eligible to reenlist 
will choose to do so. The values of R^(x^^) have been esti- 
mated from responses to previous bonus offers; more sophis- 
ticated methods for determining these functional values, for 
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example, by including the effect of past and present 
economic conditions, are currently being studied [Ref. 5]. 

As stated previously, otherwise identical deviations in 
different cells are not viewed equally. Generally it costs 
the Marine Corps considerably more to replace or train an 
individual in a highly technical MOS than in a less tech- 
nical MOS . Further, the actual number of Marines in a cell 
influences the acceptability of shortages. For example, a 
shortfall of 5 in a cell containing 1000 could be more 
easily tolerated than that same shortfall in a cell 
containing only 10. 

Therefore, each cell i is assigned the weighting coeffi- 
cient is the "training index" for cell i which 
is just the cost of training a Marine in cell i normalized 
by the maximum training cost over all cells. is the 
reciprocal of the actual number of Marines in cell i. is 
an exogenous weighting factor which permits the user of the 
model, if desired, to exercise his judgement in weighting 
cells or blocks of cells; presently, all are set to 1. 

The form of the resource function rj^(x^) is straightfor- 
ward. This function represents the expected number of 
dollars that must be allocated in order to offer cell i an 
SRB multiplier of x^ and is the product of the expected 
number of reenlistments in cell i that would result from 
offering a bonus x^, R^(x^)Ej^, and the size of the bonus 
min { ( .75)P^S^x^, $12,000} where $12,000 is 75% (the 
percentage of the bonus that the individual receives as a 
lump sum payment) of the maximum permissible bonus amount of 
$16,000. 

Problem Pi is recognizable as a nonlinear knapsack 
problem in which cells represent commodities and the budget 
represents the weight or cube constraint. However, the 
nonlinearity in the objective function and budget constraint 
make solution by standard branch and bound techniques 
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impossible [Ref. 6], and make solution by dynamic program- 
ming inefficient since standard reduction techniques 
[Ref. 7] are inapplicable. A straightforward dynamic 
programming solution is possible using the cells as stages 
and dollars as the sole state variable. However, signifi- 
cant computational difficulties would arise due to the large 
number of stages (on the order of 10'^) combined with a state 
variable that for FY86 can achieve any integer value on the 
interval [0,7x10^]. The solution methodology that is used 
could be applied directly to the nonlinear knapsack formula- 
tion. However, the explanation of the implementation of this 
methodology will be more transparent if the problem is first 
converted to a linear integer problem, in particular, a 
generalized assignment problem. The following section 
describes this new formulation. 

B. CONVERSION TO A GENERALIZED ASSIGNMENT PROBLEM 

A typical formulation for the generalized assignment 
problem is: [Ref. 8] 

Indices : 

k = 1 , . . . ,n agent 

i = 1 , . . . ,m task 



Data : 





cost if task i is assigned to agent k 
amount of resource expended if task i is 
assigned to agent k 

minimum and maximum amounts of the 
resource that may be expended by agent k 



Decision Variables: 



1 if task i is assigned to agent k 
0 otherwise 



Formulation: 



ZZ'^ik^ik 

i k 



(P2) 
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k = 1 , . . . , n 



s.t. ay. ^ ^ by. 



=1 i = 1 , . . . ,m 

k 



^ik ^ ^0>1) 

In contrast to the classical assignment problem, the 
first n constraints imply that more than one task i can be 
assigned for completion to agent k provided that the 
resource constraints on agent k are not violated. (Only one 
type of resource is allowed. ) The second set of m 
constraints requires that each task be assigned exactly 
once . 

In the SRB multiplier selection problem the model is 
simplified to the following: 



Indices : 

i=l,...,m cells 

j = 0, . . . ,n^ SRB multiplier values 



Data : 






) 


^ij 


ri(j) 


B 


budget 


^i 


maximum permissible multiplier for cell 



n^ = 5,4,3 for Zone A,B,C, respectively 
Decision Variables: 

x^j 1 if multiplier j is selected for cell i 

0 otherwise 



i 



Formulation: 



i j 



(P3a) 



s . t 



i j 



’^ij^ij 



^ B 



Z^ij = 1 

J 
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Xij e {0,1} 



Here, a single agent is required to perform all m tasks but 
has at his disposal n^ alternative processes for each. The 
above formulation is completely equivalent to PI although a 
continuous, i.e., linear programming (LP), relaxation of 
this formulation might be very weak since neither rj^(xj^) nor 
Ci(Xi) will typically be convex functions. No problems 
arise, however, since no LP relaxation is employed in the 
solution procedure. 

This formulation of the model is quite general. The 
coefficients c^j and r^j are directly calculated by evalu- 
ating the functions and rj^(x^) of the preceding 
section at all feasible values of x . Provided that the 
assumption of independence between cells is retained, alter- 
native functions are permissible. It will be seen in the 
succeeding chapter that the solution methodology is equally 
general and introduction of alternative functions is 
trivial . 
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III. SOLUTION METHODOLOGY 



With the problem and formulation now defined, this 
chapter details the approach taken to find a satisfactory 
set of SRB multipliers. Commonly, bounds for integer 
programs are established using an LP relaxation; the special 
structure of this problem and the likely weakness of the LP 
relaxation, however, suggests an alternative approach using 
Lagrangian relaxation. The relaxed formulation is extraor- 
dinarily easy to solve as an integer program and the 
resulting bound is at least as good as that produced using 
an LP relaxation. Furthermore, in the course of determining 
the bound, feasible sets of multipliers are encountered, the 
best of which is heuristically improved to yield a final 
solution. 

A. LAGRANGIAN RELAXATION TO ESTABLISH A LOWER BOUND 

The technique of Lagrangian relaxation is well suited 
for those integer programming problems that would be rela- 
tively easy to solve were it not for one or more compli- 
cating constraints. Moreover, it has been successfully and 
efficiently applied previously to generalized assignment 
problems [Refs. 8,9]. 

As described in Fisher’s excellent primer on Lagrangian 
relaxation [Ref. 10], the approach is to move complicating 
constraints into the objective function using the product of 
the Lagrangian multiplier and the constraint violation as a 
penalty term. To demonstrate, consider the linear integer 
programming problem: 

min cx (F^) 

s . t . Ax b 
Dx ^ e 

x^O and integer 
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where the constraints Ax^b are the complicating ones. This 
is relaxed to: 

min^j cx + X,(Ax-b) (P5a) 

s . t . Dx ^ e 

x2:0 and integer 



where it is assumed that this is easy to solve for a fixed 
row vector ^2:0. For such a Xt an optimal solution to 
problem P5 provides a lower bound on the objective value for 
P4. The best possible lower bound from a relaxation of the 
form of P5 is found by solving: 






min^^ 

s . t . 






cx+X(Ax-b ) 
Dx e 
X S: 0 and 






integer 



(P5b) 



Furthermore, this bound is at least as good as that provided 
by an LP relaxation [Ref. 11]. 

Examination of the generalized assignment formulation 
developed in the preceding chapter indicates that the only 
truly complicating constraint is the budget constraint; 
integer optimization with this constraint relaxed consists 
of m separable, "multiple choice" problems which are 
trivial. Thus, the relaxed formulation with scalar X is: 



max 



niirix XS^ij^ij ^^^XX^ij^ij 
i j i j 



(P3b) 



-< 



.t. ^X£j = l i = l,...,my 



Xij e {0,1} 



where the inner portion of the objective function may 
equivalently written: 



be 



minjj ^^Xij (Cij+A,rij )-A,B (eqn 3.1) 

i j 
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For a fixed Xj the inner minimization is performed by 
selecting a minimum (cj^j+Xr^^j) for each cell. The outer 
maximization problem is a convex optimization problem which 
is easily solved, since X, is a scalar variable, by first 
bracketing X, and then performing a bisection search. 

The bracket about the optimal X, is readily established. 
The lower bound is zero. For X. sufficiently large, the 
solution to the inner minimization is to expend as few 
dollars as possible, i.e., x^q = 1 for all i. , the upper 
bound on interesting values of X,s is the smallest value of X 
for which this solution is optimal. Therefore, is the 
smallest value of X satisfying: 

Cio + XriQ:^Cij+Xrij V i,jS:l (eqn 3.2) 

Recalling that r^j represents the cost of offering cell i an 
SRB multiplier j and that j=0 corresponds to no bonus, r^Q=0 
for all i. Thus, equation 3.2 becomes: 

Cio:^Cij+A.rij V i,j^l (eqn 3.3) 
and, as a result, 

L^=max{ (c^O'^^ij )/^ij ^ V i,j^l (eqn 3.4) 

For X=hjj, a tie exists between the solution Xj^Q = l for all i 
and at least one solution where x^j = l for some i and some 
j^l. For unusual data, the latter solution could be infeas- 
ible. To ensure that at least one feasible solution is 
obtained, in practice is replaced by (1 + e)Ljj where e >0. 

The solution to P3b provides a lower bound on the 
optimal solution to the SRB multiplier problem. Upper bounds 
are provided by feasible sets of SRB multipliers encountered 
while solving P3b . Such solutions arise since, for X 
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sufficiently large, the entire budget will not be consumed. 
The best of the feasible sets of SRB multipliers, that is, 
the one with the smallest objective function value and, 
correspondingly, the lowest upper bound value, is improved 
using the heuristic procedure described in the following 
section . 

B. MARGINAL RATE OF RETURN HEURISTIC 

It seems reasonable to expect that the best feasible 
solution obtained in the process of optimizing P3b would 
require the expenditure of nearly the entire SRB budget. 
There might, however, be some budget dollars remaining, a 
residual that could be apportioned to further reduce the 
objective value of the best feasible solution which is the 
present upper bound. Accordingly, a heuristic method for 
allocating this residual was developed using the concept of 
marginal rate of return. 

Given a best feasible set of SRB multipliers from the 
Lagrangian procedure, the incremental cost for each cell 
that would be realized by increasing that cell's multiplier 
by 1 is calculated for all those cells with multipliers not 
already at their maximum values. If this incremental amount 
is less than the unallocated portion of the budget it is 
possible to increase that cell's multiplier by 1 while 
remaining feasible. Let such a cell have x^j=l. Then, the 
marginal rate of return for the cell is defined to be: 

ROR^=(Cj^j-c^ (r^ ) (eqn 3.5) 

This/quotient expresses the improvement realized in the 
objective function per dollar spent when cell i's SRB multi- 
plier is increased. In the heuristic procedure, the cell 
with tne maximum positive ROR^^ has its multiplier incre- 
mented. The residual portion of the budget is reduced and 
the process is repeated until it is no longer possible to 
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increase the multiplier of any cell. The resulting set of 
SRB multipliers is accepted as the final solution and upper 
bound for the SRB multiplier problem provided that the ratio 
between the lower and upper bound is close to 1. 

Other heuristic improvement procedures are possible but 
were not implemented in this study. For example, by simul- 
taneously incrementing one cell's SRB multiplier by 1 and 
decrementing another's (an 0(m'^) operation) it might be 
possible to improve on the upper bound while remaining 
feasible. In fact, the problem could be solved to optimality 
using a branch and bound algorithm [Ref. 9] though the large 
number of decision variables might make this a slow process. 
Pursuing strict optimality, with the resultant increase in 
algorithmic complexity and solution times was not deemed 
worthwhile in light of the consistently good solutions that 
have been obtained. 
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IV. IMPLEMENTATION AND COMPUTATIONAL RESULTS 



A FORTRAN 77 computer program was written to implement 
the solution methodology described in the preceding chapter. 
Using data obtained from the Manpower Department at 
Headquarters, USMC (HQMC), the model was run on both an IBM 
3033AP and an IBM PC. To judge the robustness of the model, 
additional testing was performed using randomly perturbed 
data . 

A. COMPUTER PROGRAM 

The program that performs the Lagrangian bounding proce- 
dure and which yields the heuristically improved final solu- 
tion was written in ANSI standard FORTRAN 77. No functions 
unique to the system were used. When the program was run on 
the IBM 3033AP it was compiled by the IBM VS FORTRAN 
compiler at 0PT(2); on the IBM PC it was compiled using 
Ryan-McFarland' s RM/FORTRAN compiler with code optimization. 
The program requires six input files: Parameter file. 
Reenlistment Plan file. Reenlistment Rates file. Training 
Cost file. Cell Size file, and a Cell Weight file; File 
specifications are detailed in Appendix D. As output, the 
program produces a file containing an MOS/Zone listing with 
the proposed SRB multipliers. For comparison purposes, the 
objective value as well as the percent deviation between the 
final solution's objective value and the lower bound on the 
optimal objective are output. 

Under some circumstances the user may wish to preassign 
the SRB multipliers in specific cells. For example, the 
user may wish to declare the multipliers for all three zones 
in MOS 2112 to be some specified value. This may be accom- 
plished by "tagging” those cells in the Cell Weight file and 
entering the desired SRB multiplier in the appropriate 
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column in the data file. Additional details are provided in 
Appendix D. Comparison of the objective values before and 
after the introduction of these preassigned multipliers 
enables the user to gauge the effect of a priori multiplier 
assignments . 

B. TEST DATA 

The FY86 data that was used to develop and test the 
model was obtained from the Enlisted Plans Section of the 
Manpower Department at HQMC . Model input consists of those 
values listed in equation 2.1 and 2.2 organized into five 
data files and one problem parameter file. Because the rele- 
vant data is extracted from a variety of sources, it was 
deemed desirable to place them in different files. The 
information contained within the Parameter file is primarily 
run specific. Samples from the various input files are 
contained in Appendix B. 

C . COMPUTATIONAL RESULTS 

The data described above leads to a problem with 979 
cells. This in turn results in a problem, in the form of 
P3b, with 4895 binary variables and 980 constraints. The 
objective function value from the final solution is within 
.01% of the lower bound established through the Lagrangian 
procedure. In the process, 99.99% of the budget was 
expended. Appendix C contains a sample from the output 
f ile . 

Running time on the IBM 3033AP was 1.7 seconds. 1.4 of 
the 1.7 seconds is used to read in the data and calculate 
the coefficients Cj^j and r^j . The Lagrangian procedure and 
the heuristic use most of the remaining .3 seconds. The 
time required to write the solution was negligible. 

After preliminary development and testing on the IBM 
3033AP, the program and data were copied onto a floppy disk- 
ette and, using the RM/FORTRAN software, installed on an IBM 
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PC configured with 512Kb RAM, an 8087 Math Coprocessor and 
two 360Kb floppy disk drives. No changes to the code were 
required in moving it to the PC. On the IBM PC, the program 
runs in approximately 4 minutes. Most of this 4 minutes is 
consumed reading in the data from the floppy diskette; 
slightly faster times might be expected using a system with 
a hard disk arrangement . 

Additional testing was conducted to verify the robust- 
ness of the solution procedure and to confirm that the 
excellent results that were observed were not merely fortu- 
nate happenstance stemming from a "good" set of data. In 
this procedure, the data read by the model was perturbed in 
this manner: 

d’= Ld’"UNIF( .5,1.5 )j 
where 

d' perturbed datum 

d original datum 

UNIF(a,b) random variate generated from the continuous 

uniform distribution between a and b 

LXj largest integer not exceeding x 

This randomization process was performed on all the data 
except that response rates were capped at 1 and not discre- 
tized . 

In all, 50 additional model runs were conducted in this 
manner, each with a different data set. Solution times on 
the IBM 3033AP varied slightly, ranging between 1.5 and 1.9 
seconds. In each case, the final solution was within .02% 
of the lower bound on the optimal solution. Typically, the 
heuristic improves the best feasible solution with respect 
to the lower bound by less than .04%. In a situation in 
which exceedingly tight bounds are not necessary, the 
heuristic could, in fact, be eliminated. 



25 



V. CONCLUSIONS AND RECOMMENDATIONS 



In this thesis, the process of determining Selective 
Reenlistment Bonus multipliers subject to a limited budget 
was modeled as a mathematical program. In the course of 
doing this, an objective function based on expected devia- 
tions from cell reenlistment targets was defined. This is an 
important advance; formerly there was no method by which 
alternative sets of multipliers could be compared. Using 
this objective function, a procedure was developed which 
determines a nearly optimal set of SRB multipliers. The 
speed with which the solution is found, a modest core 
storage requirement and a compact, specialized FORTRAN code 
proves the procedure amenable for PC use. This is an impor- 
tant advantage given the PC's wide-spread availability. The 
inclusion of an exogenous weighting factor and the capa- 
bility to preset multipliers allows the user additional 
flexibility in using the model. 

The model and solution method that are employed are very 
general and only require that the assumption of inter-cell 
independence be maintained. Other penalty and resource 
utilization functions are readily incorporated into the 
model. Improvements in the estimates for the response rates 
should further increase the accuracy with which the model's 
multipliers achieve the cell targets. 

Provided the user possesses the requisite hardware (the 
IBM PC configuration explained in Chapter IV should be 
considered the minimum) and a FORTRAN 77 compiler, the 
program is completely operational. A copy of the source 
code is listed in Appendix A. Other pertinent information, 
including examples of input and output files and user 
instructions, is contained within Appendices B, C, and D. 
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The future portends increasingly restrictive budgets 
within the Department of Defense and additional pressure to 
maximize benefits given limited resources. With respect to 
maximizing the benefits realized from the SRB program, a 
nearly optimal assignment of SRB multipliers goes far 
towards achieving that goal. 
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APPENDIX A 

LISTING OF SOURCE CODE 



IMPLICIT REAL* 8 (A-H,0-Z) 

IMPLICIT INTEGER*4 (I-N) 

REAL*8 Cf 1100, 0:5 ), aTiIOO, 0:5) 

INTEGER*2 MOStllOO) .IZONEfllOO) ,IXBEST(1100) ,IX(1100) 

COMMON CINF, BUDGET, TEST, NUMCEL 

CINF=1.D22 

CALL INITAL(C , A.HLAM.MOS , IZONE , IX) 

CALL BOUNDTzL , IX , ZU . IXBEST , HLAM , C , A , COST , MOS ) 

CALL HEURIS (C , A , IXBEST , ZH , COST , MOS ) 

WRITE (16,50) 

FORMAT r lx , ' **VoVV.*******Vj^V*******V^V.-*VfVr*VrVf**Vr » ) 

PCTDEV= 100 . DO* ( 1 . DO- ZL/ ZH) 

WRITE (16,100)ZH 

FORMAT fix 'OBJECTIVE FUNCTION VALUE IS:',F11.3) 

WRITE (16, 125jPCTDEV 

FORMAT fix, 'THIS SOLUTION IS WITHIN ',F11.8, 

•' % OF THE OPTIMAL') 

WRITE (16,50) 

WRITE U6,300) 

DO 150 1=1, NUMCEL 

WRITE (16 ,200)MOS(I) , IZONE (I ) , IXBEST ( I ) 

CONTINUE 

FORMAT 7iX,I5 .3X,I2,9X,I2) 

FORMAT (3X, 'm6s' , 3X , 'ZONE ,5X, 'MULTIPLIER' ) 

STOP 

END 

SUBROUTINE INITAL ( C , A , HLAM , MOS , IZONE , IX ) 

IMPLICIT REAL* 8 (A-H,0-Z) 

IMPLICIT INTEGER*4 (I-N) 

REAL*8 A0S(3),APAY(3),C( 1100,0:5) ,A(1100, 0:5) ,R(0: 5) 
INTEGER*2 MOS U) , IZONfeTl) , IXU) 

COMMON CINF, BUDGET, TEST, NUMCEL 
READ (-10, 100)BUDGET 
10 ,200)NUMCEL 
10,100)TIMAX 

10, 

10, 

10,' 

10 ' 

HLAM=-CINF 
SCALE=1000.D0 
BUDGET = BUDGET/ SCALE 
BNSMAX=BNSMAX*PCTGFY 
TEST=CINF/1. IDO 

NUMCEL) 



READ 

READ 

READ 

READ 

READ 

READ 

READ 



300)(AOS(I),I=1,3) 
300uapayTi) ,i=i,3) 
100)BNSMAX 
400)PCTGFY 
400)Q 



READ(11,500)(MOS(I) ,IZONE(I) ,IX(I) , A (1 , 0 ) , 1= 1 , : 
READ ' ( l4 . 6 0 02 ( A ( i , 1 , A ( I , b i = 1 ^ N^ck 
DO 50 1=1, NUMCEL 

“'L=6-IZONE(I) 

(15 ,700) (C(I, J) , J=0,MAXMUL) 

,800) (A(I, 3) ,1=1, NUMCEL) 

L~ D 

E = A(1,1 



MAXMUL=6-IZONE(I 
READ 
CONTINUE 
READ (12,800 

READ U3 

DO 150 

MAXMUL=klZONE(I) 

wt=aTi,o) 

E= ■ ' 

TE ^ 



ACT=i(i,]; 
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TI=A(I,4) 

IF (TI .LT. l.DO) TI=50.D0 
DO 75 J=0,MAXMUL 
R(J)=C(I, J) 

CONTINUE 

AVGPAY=APAY(IZONE(I)) 

avgaos=aos?izoneTi ) ) 

IF (ACT .GT. O.ODO) THEN 
CW= WT ( TI / TIMAX ) / ACT 

ELSE 

CW=WT''TI /TIMAX 
END IF 

IF (E .LT. l.DO) E=1.D0 
DO 125 J=0,5 

IF Tj .le. maxmul) then 
DIFTRM=TE-E“R7j) 

IF (DIFTRM .LE. O.DO) THEN 

C ( I , J ) = CW*Q* ( DIFTRM*-'2 ) 

ELSE ^ ^ 

C ( I , J ) = CW* ( DIFTRM**2 ) 

END IF 

BONUS =PCTGFY-AVGPAY*AVGAOS-J 
IF (BONUS .GT. BNSMAX) BONUS=BNSMAX 
: . J ) = r7 J ) *E " BONU S / S GALE 
(J .GT. 0) THEN 






END IF 
ELSE 



C(I,J)=CINF 
AiI,J)=CINF 



END IF 
CONTINUE 
CONTINUE 
FORMAT (F12.0) 
FORMAT (14) 
FORMAT (3F7 .0) 
FORMAT (F4. 2) 



FORMAT (I5,2X,I1,2X,I1,2X,F4.2) 
FORMAT (8X,iFl6.0 3 




,A(110O,O: 1) 
3estIiI,mos(1) 



RETURN 

END 

SUBROUTINE BOUND ( ZL , IX , ZU , IXBEST , HLAM , C , A , COST , MOS ) 
IMPLICIT REAL*8 (A-H,0-Z) 

IMPLICIT INTEGER-4 (I-N) 

REAL-8 C( 1100, 0:1),/"”'^' 

INTEGER*2 IxTi),IXBI 
COMMON CINF, BUDGET, TEST, NUMCEL 
EPS=HLAM71.D05 
ENDL=O.DO 
ENDR= .01D0*HLAM 
ZUBEST=CINF 

AMBDA= ( ENDL + ENDR ) / 2 . DO 

CALL MINFCN (AMBDA , ZL , IX , ZU , COST , C , A , MOS ) 

IF (COST .LE. BUDGET) THEN 
ENDR= AMBDA 

IF (ZU .LE. ZUBEST) THEN 
ZUBEST=ZU 
FLAM = AMBDA 
END IF 
GOTO 100 
END IF 
ENDL=AMBDA 
ENDR=10.01D0*ENDR 
AMBDA= ( ENDL+ ENDR ) / 2 . DO 

CALL MINFCN (AMBDA , ZL , IX , ZU , COST , C , A , MOS ) 
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IF (COST .LE. BUDGET) THEN 
ENDR=AMBDA 

IF (ZU .LE. ZUBEST) THEN 
ZUBEST=ZU 
FLAM=AMBDA 
END IF 
GOTO 100 
END IF 
ENDL=AMBDA 
ENDR= 10 . 0 IDO * ENDR 
100 AMBDA= (ENDL+ENDR)/2.D0 

CALL MINFCN ( AMBDA , ZL , IX , ZU , COST , C , A , MOS ) 

IF (COST .LE. BUDGET) TH^N 
ENDR= AMBDA 

IF (ZU .LE. ZUBEST) THEN 
ZUBEST=ZU 
FLAM=AMBDA 
END IF 
FT «;F 

ENDL= AMBDA 
END IF 

IF ((ENDR-ENDL) . GT . EPS) GOTO 100 
C CALCULATE BEST FEASIBLE SOLUTION 

CALL MINFEA ( FLAM , ZU , IXBEST , COST , C , A , IX , MOS ) 

RETURN 

END 

SUBROUTINE MINFCN (AMBDA ,ZL , IX , ZU , COST , C , A , MOS ) 
IMPLICIT REALMS (A-H,0-i) 

IMPLICIT INTEGER-4 (I-N) 

REALMS CfllOO.O: 1 ) ,aTi10O,O: 1) 

INTEGER*^ IxTn,MOS(lT 

COMMON CINF, BUDGET, TEST, NUMCEL 

C0ST=0.D0 

CELTOT=O.DO 

ZU=0.D0 

DO 400 1=1, NUMCEL 

IF (MOS (I) .LT. 0) THEN 

CELMIN= C ( I , IX ( I ) ) + AMBDA*A ( I , IX ( I ) ) 

COBJ=C(i:i^(I)) 

INDEX=IX(I) 

GOTO 300 
END IF 
CELMIN=CINF 
INDEX =0 
DO 200 J=0,5 

IF (C(I,J) .LE. TEST) THEN 
Z J = C] I , J ) + AMBDA-A ( I , J ) 

CEL0BJ=C(I, J) 

IF (ZJ .LE. CELMIN) THEN 
IX(I)=J 
CELMIN=ZJ 
COBJ=CELOBJ 
INDEX=J 
END IF 
END IF 
CONTINUE 

CELTOT=CELTOT+ CELMIN 
ZU=ZU+COBJ 

COST=COST+A(I , INDEX) 

CONTINUE 

ZL = CELTOT - AMBDA*BUDGET 
RETURN 
END 

SUBROUTINE MINFEA ( FLAM , ZU , IXBEST , COST , C , A , IX , MOS ) 
IMPLICIT REALMS (A-H.O-Z) 

IMPLICIT INTEGER-4 (I-N) 

REAL-8 C( 1100,0: 1) ,aTi10O,O: 1) 

INTEGER*2 IXBEST(l) ,IX(1) ,M0S(1) 



100 



200 

300 



400 



30 



100 



200 

300 

400 



100 



200 



300 



9x( 



I) 



COMMON CINF, BUDGET, TEST, NUMCEL 
COST=O.DO 
CELTOT=O.DO 
ZU =0.D0 

DO 400 1=1, NUMCEL 

IF (MOsll) .LT. 0) THEN 
COBJ=C(I,IX(I)) 

INDEX=IX?^' 

ixbestTi 
GOTO 300 
END IF 
CELMIN=CINF 
INDEX=0 
DO 200 J=0,5 

IF (C(I,J) .LE. TEST) THEN 
ZJ = CXI , J)+FLAM''A(I , J) 

CEL0BJ=C(I, J) 

IF (ZJ .LE. CELMIN) THEN 
IXBEST(I)=J 
CELMIN=ZJ 
COBJ=CELOBJ 
INDEX=J 
END IF 
END IF 
CONTINUE 
ZU=ZU+COBJ 

COST=COST+A(I , INDEX) 

CONTINUE 

RETURN 

END 

SUBROUTINE HEURIS f C , A , IXBEST , ZH , COST , MOS ) 
IMPLICIT REAL“8 (A-H.O-Z) 

IMPLICIT INTEGER^ (1-N) 

REALMS Cf 1100,0: 1 ) , aTiiOO , 0 : 1 ) , UARRAY ( 1100 ) 
INTEGER*2 IXBEST(l) ,M0SU3 
COMMON CINF, BUDGET, TEST, NUMCEL 
COSTl=COST 
RESID=BUDGET-C0ST1 
DO 200 1=1, NUMCEL 

IF (MOS]l) .LT. 0) THEN 
UARRAY{I)=-CINF 
GOTO 2 OO 
END IF 

IF (IXBEST(I) .GE. 5 
' .GT. TEST) THEN 

UARRAY(I)=-CINF 
GOTO 2 OO 
END IF 



OR. C(I, (IXBEST(I)+1)) 




UARRAYfl)=-CINF 
GOTO 2 OO 

ELSE 

UARRAY ( I ) =RNUM/RDENOM 
END IF 

IF (UARRAY(I) .LT. 0 . DO ) UARRAY ( I )=- UARRAY ( I ) 
CONTINUE 
UBEST=-CINF 
DO 300 1=1, NUMCEL 

IF (UAR:^Y(I) .GT. UBEST) THEN 

UBEST=UARRAYU) 

INDEX=I 
END IF 
CONTINUE 

IF (UBEST .LT. (-TEST)) GOTO 400 



08) THEN 



cost 1=C0ST1-A( INDEX, IXEEST^NDEX)) 
*+A(INDEX, UXBEST(INDEX)+1) ) 
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IXBEST(INDEX)=IXBEST (INDEX) +1 
IF (COSTl .LE. BUDGET) GOTO 100 
400 ZH = O.ODO 
COST=O.DO 
DO 500 I=1,NUMCEL 



ZH=ZH+C( 

COST=COS 




500 CONTINUE 



RETURN 

END 
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APPENDIX B 



EXAMPLES OF INPUT FILES 
PARAMETER FILE 



70000000 

979 

72000 
4 5 

1010 1229 

16000 

0.75 

1.00 



5 

1340 
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REENLISTMENT RATE FILE 



1833 


1 


0.10 


0.20 


1833 


2 


0.65 


0.52 


1833 


3 


0.85 


0.90 


2111 


1 


0.24 


0.30 


2111 


2 


0.50 


0.51 


2111 


3 


0.85 


0.90 


2112 


1 


0.24 


0.30 


2112 


2 


0.50 


0.51 


2112 


3 


0.85 


0.90 


2131 


1 


0.24 


0.30 


2131 


2 


0.50 


0.51 


2131 


3 


0.85 


0.90 


2142 


1 


0.20 


0.21 


2142 


2 


0.50 


0.51 


2142 


3 


0.85 


0.90 


2144 


1 


0.20 


0.27 


2144 


2 


0.50 


0.51 


2144 


3 


0.85 


0.90 


2145 


1 


0.20 


0.21 


2145 


2 


0.50 


0.51 


2145 


3 


0.85 


0.90 


2146 


1 


0.24 


0.30 


2146 


2 


0.50 


0.51 


2146 


3 


0.85 


0.90 


2147 


1 


0.24 


0.30 


2147 


2 


0.50 


0.51 


2147 


3 


0.85 


0.90 


2161 


1 


0.20 


0.30 


2161 


2 


0.50 


0.51 


2161 


3 


0.85 


0.90 


2171 


1 


0.20 


0.30 


2171 


2 


0.50 


0.51 


2171 


3 


0.85 


0.90 


2172 


1 


0.20 


0.27 


2172 


2 


0.50 


0.51 


2172 


3 


0.85 


0.90 


2311 


1 


0.23 


.25 


2311 


2 


0.53 


0.70 


2311 


3 


0.85 


0.90 


2336 


-1 


.21 


.22 


2336 


2 


0.86 


. 900 


2336 


3 


0.85 


0.90 


2512 


1 


0.21 


0.39 


2512 


2 


0.48 


. 75 


2513 


1 


0.21 


0.39 


2513 


2 


.48 


. 75 


2519 


1 


0.21 


0.39 


2519 


2 


0.77 


. 78 


2519 


3 


0.85 


0.90 


2531 


1 


0.21 


0.39 



.21 


0 


.34 


.43 


.5 


.69 


0 


.88 


. 9 




. 90 


0 


. 90 






. 36 


0 


.63 


0.67 


. 7 


.70 


0 


.80 


0.80 




.90 


0 


.90 






. 36 


0 


.60 


0.67 


.7 


.70 


0 


.80 


0.80 




.90 


0 


. 90 






.63 


0 


. 70 


0.80 


.85 


.70 


0 


.80 


0.80 




.90 


0 


. 90 






.24 


0 


. 39 


0.53 


.63 


.70 


0 


.80 


0.80 




.90 


0 


. 90 






. 36 


0 


.48 


0.56 


0.63 


.70 


0 


.80 


0.80 




.90 


0 


. 90 






.37 


0 


.38 


0.65 


. 66 


. 70 


0 


.80 


0.80 




. 90 


0 


. 90 






. 36 


0 


.62 


0.65 


.66 


.70 


0 


.80 


0.80 




. 90 


0 


. 90 






. 36 


. 1 


630 


. 67 


. 7 


. 70 


0 


.80 


0.80 




.90 


0 


. 90 






.41 


0 


.42 


0.43 


0.63 


.70 


0 


.80 


0.80 




.90 


0 


. 90 






.31 


0 


. 32 


.43 


0.57 


.70 


0 


.80 


0.80 




. 90 


0 


. 90 






.36 


0 


.48 


0.56 


0.63 


. 70 


0 


.80 


0.80 




.90 


0 


. 90 






.26 


0 


.41 


.53 


.55 


.71 


0 


.72 


.8 




.90 


0 


. 90 






5 


0 


.94 


1.00 


1 . 


910 


.920 


. 95 




.90 


0 


.90 






.40 


0 


. 70 


. 71 


. 75 


76 


.88 


.9 




.40 


0 


. 70 


.71 


. 75 


76 


.88 


.9 




.40 


0 


. 70 


.71 


. 75 


8 


.88 


.9 




. 90 


0 


. 90 






.40 


0 


.70 


.71 


.75 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

0 

6 

6 

6 

0 
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1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 



REENLISTMENT PLAN FILE 



1 

5 

9 

31 

2 

4 

4 

2 

3 
42 

4 
9 

10 

3 
1 

5 
1 
1 
8 
1 
1 
9 
1 
1 

10 

1 

2 

41 

30 

2 

20 

4 
1 

34 

6 

3 
49 

4 
3 
1 
1 
3 

3 

4 
6 
1 
1 
1 
6 



22 

29 

12 

0 

8 

15 
10 
20 
13 
94 

0 

6 

0 

2 

6 

0 

0 

0 

0 

0 

4 
7 
0 

11 

16 
6 

5 

116 

0 

4 
27 

5 
0 
0 

6 
2 
0 
0 
0 

13 

13 

1 

10 

21 

9 

2 

9 

3 

6 
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0121 

0121 

0121 

0131 

0131 

0131 

0151 

0151 

0151 

0161 

0161 

0161 

0193 

0193 

0193 

0211 

0211 

0211 

0231 

0231 

0231 

0241 

0241 

0241 

0251 

0251 

0251 

0311 

0311 

0311 

0313 

0313 

0313 

0331 

0331 

0331 

0341 

0341 

0341 

0351 

0351 

0351 

0352 

0352 

0352 

0369 

0369 

0411 

0411 



CELL SIZE FILE 



1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 



3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

2 

3 

1 

2 



1098 
292 

4 

1151 
301 
11 
3567 
1036 
102 
194 
104 
38 
14 
422 
640 
18 
57 
24 
256 
221 
59 
24 
66 
30 
24 
76 
21 

1 13370 

2 1010 
34 

252 
62 
2 

2344 
252 

7 

2600 
238 
6 

2420 
245 

8 

1034 
252 

5 

1010 
831 
419 
176 
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TRAINING COST FILE 



0844 


2 


5211 


0844 


3 


5211 


0847 


1 


14027 


0847 


2 


14027 


0847 


3 


14027 


0848 


2 


15600 


0848 


3 


15600 


0861 


1 


16000 


0861 


2 


16000 


0861 


3 


16000 


1141 


1 


3256 


1141 


2 


3256 


1141 


3 


3256 


1142 


1 


9865 


1142 


2 


9865 


1142 


3 


9865 


1161 


1 


3484 


1161 


2 


3484 


1161 


3 


3484 


1169 


3 


12775 


1171 


1 


4386 


1171 


2 


4386 


1171 


3 


4386 


1181 


1 


3500 


1181 


2 


3500 


1181 


3 


3500 


1182 


1 


3500 


1182 


2 


3500 


1182 


3 


3500 


1316 


1 


7041 


1316 


2 


7041 


1316 


3 


7041 


1341 


1 


4536 


1341 


2 


4536 


1341 


3 


4536 


1345 


1 


8853 


1345 


2 


8853 


1345 


3 


8853 


1349 


3 


25216 


1371 


1 


2826 


1371 


2 


2826 


1371 


3 


2826 


1391 


1 


6705 


1391 


2 


6705 


1391 


3 


6705 


1411 


1 


4386 


1411 


2 


4386 


1411 


3 


4386 


1431 


1 


5311 
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2841 

2841 

2841 

2861 

2861 

2861 

2871 

2871 

2871 

2874 

2874 

2874 

2875 

2875 

2875 

2881 

2881 

2881 

2882 

2882 

2882 

2884 

2884 

2884 

2885 

2885 

2885 

2886 

2886 

2886 

2887 

2887 

2887 

2889 

2889 

2889 

3043 

3043 

3043 

3044 

3044 

3044 

3051 

3051 

3051 

3052 

3052 

3052 

3061 



CELL WEIGHT FILE 



1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
2 0 1.00 
3 0 1.00 

1 0 1.00 
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APPENDIX C 



SAMPLE FROM OUTPUT FILE 



Vf Vr Vr Vr Vc Vr Vc* Vf Vr Vr Vc Vr Vr Vr Vr Vc V»- ' V Vf Vr Vr Vc Vr 'V ^ 

OBJECTIVE FUNCTION VALUE IS: 1768.541 

THIS, SOLUTION. IS. WITHIN, 0.00006802 % OF THE OPTIMAL 

Vr^r^V^.'V,''.V'!V';ViV';j';V?V':V'5V':V'5r'5V':V'5V'5V'5V*'!V'>V'',V'5V';V',V','.'':V':VVc’ 

MOS ZONE MULTIPLIER 

121 1 1 

121 2 0 

121 3 1 

131 1 1 

131 2 0 

131 3 1 

151 1 2 

151 2 0 

151 3 1 

161 1 5 

161 2 0 

161 3 0 

193 1 5 

193 2 4 

193 3 0 

211 1 5 

211 2 4 

211 3 1 

231 1 5 

231 2 0 

231 3 1 

241 1 5 

241 2 0 

241 3 0 

251 1 5 

251 2 4 

251 3 1 

311 1 1 

311 2 4 

311 3 1 

313 1 5 

313 2 0 

313 3 1 

331 1 5 

331 2 4 

331 3 1 

341 1 1 

341 2 0 

341 3 1 

351 1 0 

351 2 0 

351 3 0 

352 1 0 

352 2 0 

352 3 1 

369 2 4 

369 3 0 
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APPENDIX D 



USER INSTRUCTIONS 



A. DATA FILES 

Input to the program is provided through six data files. 
File names , specif ications and entries are detailed in the 
accompanying tables. In general, data values should be 
right justified within the fields. Rates, percentages, and 
the exogenous weighting factors may be placed anywhere 
within the specified field provided they are entered as two 
place decimals (e.g., .25, .03 or 1.15). With the exception 
of the Parameter file, all files contain the same entries in 
the first 8 columns. The program reads data values sequen- 
tially from the files. Therefore, it is vital that the 
MOSs , Zones and accompanying entries from one file corre- 
spond line for line with all other files. For example, if 
line 15 from one file contains entries pertaining to MOS 
0193, Zone C, then line 15 from all other files should refer 
to MOS 0193, Zone C. 

Those cells for which ‘ a non- zero SRB multiplier is 
expressly prohibited, such as the cells in the Band and 
Marine Corps Exchange MOSs , should be eliminated from the 
data files. Similarly, there are MOSs held only by very 
senior Marines (E-8's and E-9's) not eligible for bonuses 
due to years-of-service requirements exceeding 14 years. 
These, too, should be purged from the data files. As 
written, SRBMULT . FOR can read a maximum of 1100 cells. If at 
some point it becomes necessary to increase this, it will 
require changing the array and matrix dimensions within the 
program from 1100 to the desired number. 



40 





TABLE la 

DATA FILE SPECIFICATIONS 
Reenlistment Plan File (RPLAN.DAT) 


Line 


Column 


Data Description 


all 


1-5 

6-7 

8 

9-18 

19-28 


MOS 

blank 

Zone (1 for A, 2 for B, 3 for C) 
Number of Marines eligible for 
reenlistment in cell 
Reenlistment target for cell 







TABLE lb 






DATA FILE SPECIFICATIONS 






Cell 


Size File (ACTNUM.DAT) 




Line 


Column 


Data Description 




all 


1-5 

6-7 

8 

9-14 


MOS 

blank 

Zone (1 for A, 2 for B, 
Total number of Marines 


3 for C) 
in cell 
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TABLE Ic 

DATA FILE SPECIFICATIONS 
Training Cost File (TCOST.DAT) 


Line 


Column Data Description 


all 


1-5 MOS 

6-7 blank 

8 Zone (1 for A, 2 for B, 3 for C) 

9-19 Training cost (dollars) 





TABLE Id 

DATA FILE SPECIFICATIONS 
Response Rate File (RRATE.DAT) 


Line 


Column Data Description 


all 


1-5 MOS 

6-7 blank 

8 Zone (1 for A, 2 for B, 3 for C) 

9-14 Response to SRB multiplier of 0 

15-20 Response to SRB multiplier of 1 

21-26 Response to SRB multiplier of 2 

27-32 Response to SRB multiplier of 3 

33-38 Response to SRB multiplier of 4 

(if applicable) 

39-44 Response to SRB multiplier of 5 

(if applicable) 
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TABLE le 

DATA FILE SPECIFICATIONS 
Cell Weight File (WEIGHT.DAT) 



Line Column Data Description 



all 



1-5 MOS 

6-7 blank 

8 Zone (1 for A, 2 for B, 3 for C) 

9-10 blank 

11 Preset SRB multiplier (0 if not 

E reset ) 
lank 

14-18 Exogenous weighting factor (typic- 
ally, use 1.00. If other values are 
used, include decimal, e.g. .50) 



TABLE If 

DATA FILE SPECIFICATIONS 
Parameter File (PARAM.DAT) 



Line Column Data Description 



1 

2 

3 

4 



5 



6 

7 

8 



1-12 

1-4 

1-12 

1-7 

8-14 

15-21 

1-7 

8-14 

15-21 

1-12 

1-4 

1-4 



Budget available for current year 
bonuses (dollars) 

Total number of cells 

Maximum training cost (dollars) 

Zone A avg reenlistment duration 
Zone B avg reenlistment duration 
Zone C avg reenlistment duration 
(years ) 

Zone A avg monthly base pay 
Zone B avg monthly base pay 
Zone C avg monthly base pay 
(dollars J 

Maximum permissible bonus (dollars) 
Percentage of bonus awarded as lump 
sum (in decimal form, e.g., .75) 

Over/under factor 
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B. RUNNING THE PROGRAM 



To use the program, the user should possess, at a 
minimum, the IBM PC configuration detailed in Chapter IV and 
a FORTRAN 77 compiler. The program was developed and tested 
using the Ryan-McFarland RM/FORTRAN compiler. However, the 
code used in the program is sufficiently generic that other 
FORTRAN 77 compilers should be acceptable. 

A diskette ("program diskette") containing the following 
files will be provided to the Enlisted Plans Section, 
Manpower Department, HQMC : HOWTO.DOC, RUNCHK.BAT, 
FILCHK.FOR, FILCHK.EXE, RUNSRB.BAT, SRBMULT.FOR, and 
SRBMULT.EXE. The file HOWTO.DOC contains a copy of oper- 
ating instructions. Files with filetype .BAT are executive 
routines, those with filetype .FOR are FORTRAN source code 
and files with filetype .EXE are executable programs. Use of 
these programs is detailed in the following sections. Prior 
to running any of these programs certain parameters within 
the DOS configuration file CONFIG.SYS must be changed. If 
this file does not already exist, it will be necessary to 
create it. The RM/FORTRAN that was used to compile 
SRBMULT.EXE and FILCHK.EXE requires that the number of files 
that can be opened concurrently be increased to at least 
ten. Therefore, the following line must be added to 
CONFIG.SYS: 



FILES = 10 

Similarly, the number of disk buffers allocated by DOS at 
system startup must be increased to ten. Accordingly, add: 

BUFFERS = 10 

to CONFIG.SYS. Once these changes have been made to 
CONFIG.SYS, reboot DOS. 
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1 . Checking Data Files for Conformity 

The FORTRAN program FILCHK.EXE is provided to check 
the files for MOS and Zone conformity prior to using 
SRBMULT.FOR. To use FILCHK.EXE on a system with two disk 
drives, place the program diskette in B drive and the disk- 
ette containing the data files ("data diskette") in A drive. 
Enter the following DOS commands: 

Path = B:\ 

A: 

The first command will cause DOS to check the B drive for 
executable files. The second sets the A drive as the default 
drive. At the system prompt, enter RUNCHK to execute the 
program RUNCHK.BAT. This executive routine will set neces- 
sary file definitions and then execute FILCHK.EXE. On a 
system with a hard disk, place all files in the same direc- 
tory on the C drive (hard disk) and then execute RUNCHK.BAT 
as before. 

The FILCHK program uses the sequence of MOSs and 
Zones in the file WEIGHT.DAT as a reference set and checks 
the other four data files for conformance to it. If a non- 
conforming entry is discovered,' the applicable file defini- 
tion number and line number are written to the screen. The 
remaining files are then checked. Once detected errors are 
corrected, RUNCHK.BAT should be executed again until all 
files conform with the reference as evidenced by no error 
returns . 

2 . Running SRBMULT 

To run SRBMULT.EXE on a system with two disk drives, 
place the program diskette in the B drive and the data disk- 
ette in the A drive. Enter the following DOS commands: 

PATH = B:\ 

A: 

RUNSRB 
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Several minutes later, an output file MULPLN.DAT will be 
created on the diskette in the A drive. It will contain an 
MOS and Zone listing, the recommended SRB multipliers, the 
objective function value and the percent deviation from the 
optimal. The objective function value is useful in judging 
the effect of introducing exogenous weighting factors or 
presetting multipliers. Instructions on how to do both of 
these is contained in the following section. 

Running the program is somewhat simpler on a system 
with a hard disk. All files should be placed in a single 
directory on the C drive (hard disk) and that drive should 
be made the default drive. SRBMULT.EXE is initiated by 
typing RUNSRB to execute RUNSRB.BAT. 

C. EXOGENOUS WEIGHTING FACTORS, PRESET MULTIPLIERS AND THE 

OVER/UNDER FACTOR 

Inclusion of an exogenous weighting factor (columns 
14-18 in WEIGHT.DAT) permits the user to exercise additional 
control in weighting the effects of reenlistment shortfalls 
in particular cells based on his judgement. Weighting 
factors Wj^ greater than the default value of 1.0 result in 
comparatively greater penalties for failing to meet reen- 
listment targets. Conversely, weighting factors less than 
1.0 result in lesser penalties. The effect of the factor is 
linear. Thus a 2.0 results in a penalty twice as large as a 
1.0 and a factor of .5 reduces the penalty by one-half. As 
a weighting scheme based on training costs and manning 
levels is already utilized, this capabililty should be used 
with caution. In general, a value of 1.0 should be used for 
all cells. 

The user can also preset SRB multipliers. This might be 
necessary if factors other than those considered by the 
model dictate that particular cells be assigned specific 
multipliers. This is accomplished by entering a minus sign 
(-) in column 1 of the pertinent cell's data entry in the 
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file WEIGHT.DAT. The desired multiplier is then entered in 
column 11 of the same line. In the output file, the minus 
sign will appear next to the appropriate MOS/Zone thereby 
"flagging" those cells with predetermined multipliers. 

The last entry in PARAM.DAT is the over/under factor Q, 
the weighting factor employed to permit the user to estab- 
lish the relative importance of overshooting or under- 
shooting reenlistment targets. For example, a over/under 
factor of .70 implies that the penalty incurred in the model 
for overshooting the reenlistment target by some number is 
only 70% of the penalty that would be incurred for under- 
shooting the target by the same number. Since an overage in 
a particular cell is generally considered more acceptable 
than a shortage, the factor will typically be less than 1. 
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